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Probability of US heat waves affected by a 
subseasonal planetary wave pattern 

Haiyan Teng 1 *, Grant Branstator 1 , Hailan Wang 2 , Gerald A. Meehl 1 and Warren M. Washington 1 


Heat waves are thought to result from subseasonal atmospheric variability. Atmospheric phenomena driven by tropical 
convection, such as the Asian monsoon, have been considered potential sources of predictability on subseasonal timescales. 
Mid-latitude atmospheric dynamics have been considered too chaotic to allow significant prediction skill of lead times beyond 
the typical 10-day range of weather forecasts. Here we use a 12,000-year integration of an atmospheric general circulation 
model to identify a pattern of subseasonal atmospheric variability that can help improve forecast skill for heat waves in the 
United States. We find that heat waves tend to be preceded by 15-20 days by a pattern of anomalous atmospheric planetary 
waves with a wavenumber of 5. This circulation pattern can arise as a result of internal atmospheric dynamics and is not 
necessarily linked to tropical heating. We conclude that some mid-latitude circulation anomalies that increase the probability 
of heat waves are predictable beyond the typical weather forecast range. 


T he increasing severity of heat waves in recent decades and 
projections of even more intense and frequent heat waves in 
the future 1-3 has resulted in a growing demand for skillful 
predictions of these high-impact events beyond the 10-day forecast 
range. On the subseasonal timescale, low-frequency intraseasonal 
tropical phenomena such as the Madden-Julian oscillation 4,5 , 
fluctuations of Asian monsoon precipitation and slowly varying 
boundary conditions are commonly regarded as the primary 
potential sources of predictability. The mid-latitude circulation, 
however, is dominated by weather noise and its evolution is strongly 
chaotic making it unpredictable after about 10 days for typical 
initial conditions. However, this 10-day limit of prediction skill 
does not exclude the possibility that some circulation states can 
be substantially more predictable than the average 6 ' 7 , either from 
having large amplitude or being associated with patterns that are 
intrinsically of low frequency. 

Some case studies 8-10 have suggested that US heat waves 
(and associated droughts) may be preceded and accompanied 
by quasi-stationary large-scale mid-latitude atmospheric Rossby 
waves. However, the scarcity of extreme events in the observational 
records together with the inherent noisiness of mid-latitude weather 
fluctuations make it difficult to test this possibility. To avoid this 
problem here we use a 12,000-year simulation from an atmospheric 
general circulation model (GCM), called Community Atmospheric 
Model version 3 (CAM3), so that there are sufficient samples to 
examine these unusual events. In this simulation monthly varying 
sea surface temperature (SST) is prescribed to have present-day 
climatological values. As a result, atmospheric fluctuations are 
mainly produced by internal atmospheric dynamics. 

US heat wave statistics 

We devote our attention to events in which surface air temperature 
(SAT) is persistently high somewhere over the continental United 
States. Specifically, a day during June-August is considered to be 
a heat wave day if on that day and each of four succeeding days 
there are more than ten grid points (corresponding to 5% of the 
searched domain) over the US continental area within the domain 
of 125° W-70° W, 25° N-50° N with daily averaged SAT exceeding 


a threshold value, and the centre of these warm points does not 
move faster than 5° latitude or longitude per day. The temperature 
threshold we use varies with grid point as well as day of the year. It 
corresponds to the 97.5 percentile for historical daily temperatures 
within a 15-day window centred on the day of the year of the 
potential heat wave day. Consecutive heat wave days are grouped 
into a single heat wave event. 

We also apply the same procedure to the National Centers 
for Atmospheric Prediction/National Center for Atmospheric 
Research reanalysis 11 (NCEP/NCAR reanalysis) during the period 
of 1948-2012. A comprehensive comparison of the model and the 
observations is difficult owing to the brevity of the observations. 
Nevertheless, there are qualitative similarities between the model 
and the reanalysis data regarding the life span, intensity and 
frequency of the heat waves (Fig. 1), which indicate that the basic 
characteristics of the simulated heat waves are realistic. 

A zonal wavenumber-5 precursor pattern 

To look at the temporal evolution of planetary waves associated 
with the heat waves in the model, we designate the first day of 
a heat wave event as dayO and construct composites on each 
day from day -20 to day 5 using daily subseasonal anomalies 
of 300 hPa streamfunction (Fig. 2). Subseasonal anomalies are 
daily departures from both the long-term mean for each day 
and the seasonal (June-August) mean for any particular year. 
To avoid contamination from previous heat wave events, we 
use only 2,300 events that have no heat wave days in the 
preceding 20 days. From day -20 to day 5 (Fig. 2), there is 
a wavenumber-5 structure slowly propagating westwards as 
highlighted by the yellow dashed lines, at a speed of roughly 2-3° 
longitude per day. Both the spatial structure and the movement 
are reminiscent of atmospheric Rossby waves trapped in the mid- 
latitude jet stream waveguide 12,13 . The wavenumber-5 structure is 
equivalent barotropic, though the lower tropospheric circulation 
anomalies (figure not shown) over the Asian summer monsoon 
region are fairly weak. 

On dayO, the entire continental United States and Mexico are 
dominated by anticyclonic circulation aloft and much reduced 
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Figure 1 | Comparison of the model heat wave statistics with observations. Climatological SAT anomalies (°C; top panels) and frequency of occurrence 
(middle panels) of the US heat wave days and life span of the events (bottom panels) in the NCEP/NCAR reanalysis (left panels) and CAM3 simulation 
(right panels). 


precipitation, which is consistent with observations of heat 
waves in nature 3 ' 14,15 . Thereafter, the two positive streamfunc- 
tion centres over the Bering Strait and the continental US dis- 
appear on day 5 and day 15 (not shown) respectively and the 
wavenumber-5 pattern fades away. 

The Plumb flux 16,17 has been frequently used as a diagnostic 
of the energy propagation of stationary atmospheric Rossby 
waves. It indicates (arrows in Fig. 2) that although the phase 
of the wavenumber-5 Rossby wave propagates slightly westwards 
from day -15 to day 0, energy moves consistently eastwards 
and completes a global circle from day -10 to day 10, which is 
roughly consistent with the group velocity of a stationary zonal 
wavenumber-5 disturbance trapped in the tropospheric jet 18 . Over 
the Pacific, there are also substantial similarities in both the Plumb 
fluxes and the Rossby wave with some observed heat wave events 
(for example, the 1980 heat wave/drought 10 ) on about day 0. 

The strongest fluxes are located at the jet entrance and exit 
regions over both the Pacific and the North Atlantic oceans, 
suggesting that instabilities of the jet streams or interaction with 
transient eddies 19 may be instrumental during the life cycle of 
the heat waves. That synoptic eddies may be contributing to 
maintenance is supported by a plot (Supplementary Fig. 3) of 
the streamfunction tendency from synoptic eddy vorticity fluxes. 
That plot shows that the anomalous tendencies have the same 
sign as the streamfunction anomalies during day -20 to day 5 
in the two jet regions. 

Although there is consistent change in precipitation (shading 
in Fig. 2) as the energy propagates eastwards from day -10 to 
day 5, there are no organized precipitation anomalies in either 
the tropics or in the Asian summer monsoon region around 


day -15 or earlier. Although SST is fixed in the experiment, 
the subseasonal variability of precipitation is comparable to that 
of the observations (Supplementary Fig. 1). Therefore, we do 
not think this lack of connection with the tropics is caused by 
any deficiency of tropical precipitation anomalies resulting from 
the simple treatment of SSTs. Instead, it further supports the 
interpretation that the wavenumber- 5 structure is produced by 
internal dynamics in the mid-latitudes, with submonthly transient 
vorticity flux anomalies playing an important role, similar to the 
findings in other studies 8 . 

Pattern of intrinsic subseasonal variability 

Unlike in the model simulation there are no statistically significant 
subseasonal precursor patterns in the reanalysis fields. This is 
probably because there are so few heat waves in the short 
observational record that the signal cannot be distinguished from 
noise with statistical certainty. However, a similar wavenumber-5 
pattern has been noted in nature in boreal winter 20 and summer 21 ' 22 . 
It stands out as first empirical orthogonal function (EOF1) of 
the subseasonal variability of 300 hPa streamfunction monthly 
anomalies in the reanalysis data (Fig. 3, top). As centres of action in 
EOF patterns do not necessarily co-vary, we calculate the coherency 
in daily reanalysis streamfunction time series from the five centres. 
This calculation shows significant coherency between any pair for 
periods longer than 20 days (figure not shown), an indication that 
the summer wavenumber-5 pattern is likely to be a physical mode 
on subseasonal timescales. 

In CAM3, the wavenumber-5 pattern is represented by EOF2, 
whereas EOF1 emphasizes a zonally symmetric component of vari- 
ations in the circulation (Fig. 3, second and third panels). Pattern 


2 NATURE GEOSCIENCE | ADVANCE ONLINE PUBLICATION | www.rature.com/naturegeoscierce 

© 2013 Macmillan Publishers Limited. All rights reserved. 








1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 


NATURE GEOSCIENCE DOI: 10.1038/NGE01988 


ARTICLES 


o 

CN 

> 

ro 

Q 


un 

> 


a 


o 

> 


Q 


LD 

> 

ro 

Q 


60° N 
30° N 
0 ° 
60° N 
30° N 
0 ° 
60° N 
30° N 
0 ° 
60° N 
30° N 
0 ° 
60° N 


o 

Q 


30° N 
0 ° 


60° N 


LD 

> 

ro 

Q 


30° N 
0 ° 



0.5 

0.4 


0.3 


P 0.2 


— 0.1 


0.1 


3 

3 

Q. 


0.2 


-0.3 


-0.4 


-0.5 


Figure 2 | Temporal evolution of the heat waves. Contours represent 300 hPa streamfunction anomalies at ±0.2, ±1.0 x 10 6 m 2 s -1 levels (most areas 
with absolute values larger than 0.2 x 10 6 m 2 s -1 are significant at the 95% level; see Supplementary Figs 4, 5 for significance tests). Shading represents 
composite precipitation (mm d -1 ), with stippling indicating the 95% significance level from Student's t-test, and arrows are the Plumb flux vectors with 
magnitudes larger than 0.1 m 2 s -2 . Dashed yellow lines in the top four panels highlight the position of the five mid-latitude high-pressure areas (the 


wavenumber-5 pattern) and their movement as a function of time. 

correlations between EOF2 and the heat wave streamfunction 
composites (Fig. 2) are significant at the 99% level on each day 
from day -20 to day -12 based on a Monte Carlo test. EOF1 in 
the reanalysis seems to combine the characteristics of the leading 
two EOFs in CAM3. 

The CAM3 EOF2, with its distinctive wavenumber-5 pattern, 
has a stronger connection with US heat waves than does EOF1. 
This is seen by two-dimensional probability distribution function 
(PDF) plots of SAT averaged over the continental United States 
within 125° W-70° W, 25° N-50° N versus projections of daily 300 
hPa streamfunction anomalies onto EOF1 and EOF2 (Fig. 3 bottom 
middle and right panels). There is a clear linear relationship between 


the second principal component (PC2) and SAT when all days are 
considered. But for a given value of PC2, SAT is much stronger 
on heat wave days than it normally is, suggesting that factors 
other than PC2 must also contribute to the formation of the heat 
waves. Also presented in the bottom row of Fig. 3 is a scatter 
plot for reanalysis domain-averaged SAT and projections of daily 
300 hPa streamfuction anomalies onto the zonally asymmetric 
component of reanalysis EOF1. In this diagram the connection 
between wavenumber- 5 and SAT, including extreme events, is again 
present though not as clearly as in the large GCM data set. 

Based on a strong connection between the strength of the 
Indian summer monsoon and the wavenumber- 5 pattern in nature, 
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Figure 3 | Prominent subseasonal variability patterns, a, EOF1 of June-August subseasonal 300 hPa streamfunction monthly anomalies in the reanalysis 
(top panel) and the leading two EOFs (middle and bottom) in CAM3. b, Shown is a scatter plot of daily US continental SAT versus projections of daily 
streamfunction subseasonal anomalies on the zonally asymmetric component of EOF1 in the reanalysis (left) and probability density functions (PDFs) of 
daily US continental SAT and streamfunction projections onto the two leading EOFs in CAM3 (middle and right) for all summer days (blue) and the heat 
wave days (red). 


it has been hypothesized that an interaction with the Indian 
summer monsoon heat source is instrumental in maintaining 
the wavenumber-5 pattern 21 . For CAM3 we have examined PDFs 
of daily precipitation anomalies in both the Bay of Bengal and 
the western North Pacific, which are the key regions for the 
Asian summer monsoon. We find no distinct shifts in these PDFs 
during days when the wavenumber-5 pattern is either developing 
or at maturity. This suggests the monsoon precipitation is not 
key to producing the summer wavenumber-5 pattern in CAM3. 
Given the central role this pattern seems to play in initiating 
and maintaining US heat waves, this result is also consistent 
with our finding that tropical precipitation anomalies are not 
connected with heat wave events in this model. However, we 
cannot exclude the possibility that monsoon rainfall may initiate 
or enhance heat waves in fully coupled (atmosphere-ocean) 
experiments or in nature by stimulating or strengthening the 
wavenumber-5 structure. 


Implication for subseasonal prediction of US heat waves 

The suggestion in Fig. 2 that wavenumber-5 anomalies tend to 
precede US heat waves by 15 days makes them potentially valuable 
for subseasonal prediction of these extreme events. To investigate 
this possibility further, we project seven-day running mean 
streamfunction anomalies onto the day -15 composite pattern 
shown in Fig. 2 (the precursor pattern) and calculate the fraction 
of cases that develop into heat waves 15 days later as a function 
of the strength of the projection. The conditional probability (red 
bars in Fig. 4) indicates the stronger the amplitude of a positive 
episode of the precursor pattern the greater the likelihood of a heat 
wave two weeks later. Without any precondition the probability 
of a randomly picked day being a heat wave day is about 1.5% 
(the number of heat wave days divided by the total number of 
summer days in the simulation, denoted as P 0 ). The figure indicates 
a one standard deviation projection onto the composite pattern 
and doubles the chances of a heat wave 15 days later; a two 
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Standard deviation 

Figure 4 | Probability of the heat waves in CAM3 conditioned on the 
standardized amplitude of the wavenumber-5 pattern 15 days earlier 
(red bars). The probability of the heat waves without any precondition (Po) 
is about 0.015 (the lowest green dashed line). Asterisks and black vertical 
lines represent the average and 5th to 95th percentile range, respectively, 
of the conditional probability of 1,000 random patterns of subseasonal 
variability. The black and blue dashed lines are the PDFs of amplitude of the 
wavenumber-5 pattern in all summer days and on day -15 of the heat wave 
days, respectively. 

standard deviation projection quadruples the chances compared 
with a randomly picked case. However, initial anomalies of -0.5 
standard deviation reduce the chance of having a heat wave 15 
days later by half. A Monte-Carlo-based significance test confirms 
that such increases in conditional probability are unlikely to have 
happened by chance sampling. 

Figure 4 suggests that the wavenumber-5 pattern, through 
its influence on the likelihood of US heat waves, may benefit 
probability forecasts of these extreme events on subseasonal 
timescales. This conclusion is based on the behaviour of a GCM, 
but its applicability to nature is bolstered by the fact that the 
wavenumber-5 pattern is also a prominent subseasonal variability 
pattern in nature. Thus understanding the origins of this pattern 
is important. In particular, it needs to be recognized that this 
pattern can be generated solely by internal dynamics rather than 
being a response to Asian monsoon rainfall or SSTs. Hence 
accurate predictions of the monsoon or SSTs may not guarantee 
subseasonal predictions of US heat waves. However, our study does 
not preclude the possibility, as the work of others has suggested, 
that heat waves can also be stimulated by SST anomalies 23 or 
land surface conditions 24 . 

Methods 

Model and experiment. The 12,000-year atmospheric GCM simulation 
was conducted with the CAM3 coupled with the Community Land Model, 
which are the atmosphere and land components of the comprehensive 
atmosphere-ocean-land-sea-ice fully coupled climate model known as the 
Community Climate System Model version 3 (ref. 25). We ran CAM3 at a T42 
horizontal resolution, which corresponds to approximately a 2.8° latitude/longitude 
grid spacing. Monthly varying SST and sea-ice extent, as well as external forcing 
(solar, anthropogenic), were specified to their present-day climatological values. 
Although planetary waves produced by CAM3 have a bias in their relationship 
between frequency and zonal wavenumber (Supplementary Fig. 2), we find in 
results not described here a similar wavenumber-5 pattern leading the heat waves 
in a simulation produced by the latest, fully coupled version of the model 26 , which 
simulates a more realistic wavenumber-frequency relationship. 


Measure of the heat wave statistics. Based on our definition, there are 16,199 
heat wave days from 5,949 events in the 12,000-year simulation and 27 heat wave 
days from 17 events in the reanalysis during the period of 1948-2012. Heat wave 
intensity is represented by a composite of SAT anomalies on all heat wave days. To 
reveal the locations of the heat waves, we define heat wave frequency at a grid point 
as the number of occurrences of SAT exceeding the threshold value on a heat wave 
day divided by the total number of heat wave days. The life span of a heat wave 
event equals of the number of consecutive heat wave days plus 4. 

Significance test. Two significance tests were conducted for the wavenumber-5 
300 hPa streamfunction composite precursor patterns in Fig. 2. One is the 
standard one-sample Student f-test (Supplementary Fig. 4) and the other is a field 
significance test (Supplementary Fig. 5). For the latter, we employed a Monte Carlo 
technique in which we calculated averages of randomly drawn daily maps using the 
same sample size as in Fig. 2. Of the Fig. 2 composites the day -20 composite has the 
smallest area that passed the f-test. Of the 10,000 randomly generated composites 
only four had larger area that passed the test than the day -20 composite. Hence the 
fields in Fig. 2 composites are highly significant. 

To test the significance of the prediction results in Fig. 4 we have also 
conducted a Monte-Carlo-based test. For this test the null hypothesis is that 
increasing conditional probability of heat waves with respect to pattern amplitude 
(Fig. 4) can also be found for patterns with arbitrary structure. One thousand 
random patterns are generated from random combinations of 30 EOFs that are 
derived from 300 hPa streamfunction subseasonal variability from which variability 
associated with the wavenumber- 5 precursor pattern has been removed. As shown 
in Fig. 4, we find that the probability of heat waves is much more strongly associated 
with the amplitude of the precursor pattern than with amplitudes of these random 
patterns. We therefore reject the null hypothesis. 

Results in Fig. 4 remain valid if we apply a cross-validation method, namely, 
we use half of the 12,000-year data to construct the heat wave composite, then 
apply it to the other half of the data for finding the likelihood of a heat wave as a 
function of the amplitude of this pattern. 
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